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Abstract. We solve an extremal problem that arises in the study of the refractive indices 
of passive metamaterials. The problem concerns Hermitian functions in H 2 of the upper 
half-plane, i.e., H 2 functions satisfying f(—x) = f(x). An additional requirement is that 
the imaginary part of / be nonnegative for nonnegative arguments. We parameterize the 
class of such functions whose real part is constant on an interval, and solve the problem 
of minimizing the imaginary part on the interval on which the function's real part takes 
a given constant value. 



We consider in this note an extremal problem that arose in investigations of certain 
electromagnetic parameters of artificial materials (metamaterials). The physical interpre- 
tation of our solution in terms of bounds for refractive indices and theoretical limitations 
for the design of metamaterials is described elsewhere jHj; the purpose of the present work 
is to give an account of the underlying mathematical problem, which seems to be of some 
independent interest. 

We will be dealing with Hardy spaces H p of the upper half-plane {z = x + iy : y > 0}. 
We are primarily interested in H 2 , but it is convenient to have at our disposal the whole 
range of spaces corresponding to < p < oo. For < p < oo, H p consists of those analytic 
functions / in the upper half-plane for which 



H°° is the space of bounded analytic functions. A function / in H p has a nontangential 
boundary limit at almost every point of the real axis, and the corresponding limit function, 
also denoted /, is in LP = L P (R). Indeed, the L p norm of the boundary limit function 
coincides with the H p norm introduced above. Thus we may view H p as a subspace of L p . 
We refer to pQ for these and other basic facts about H p , as well as the twin theory of H p 
of the unit disk. (We will make a reference to the disk setting at one point.) 

The Hilbert space H 2 is the image of L 2 (M + ) under the Fourier transform. In practice, 
it is quite common that one considers functions in H 2 that are Fourier transforms of real- 
valued functions in L 2 . This leads to the following symmetry condition: f(—x) = f(x). 
Functions / satisfying this condition will be referred to as Hermitian functions. Thus 
Hermitian functions have even real parts and odd imaginary parts. 
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The Hilbert transform of a function u in LP (1 < p < oo) is defined as 




It acts boundedly on LP for 1 < p < oo and iso metrically on L 2 . If u is a real- valued 
function in IP for 1 < p < oo, then u + iu is in and so the role of the Hilbert transform 
is to link the real and imaginary parts of functions in H p . We will only work with Hermitian 
functions, and we will be interested in computing real parts from imaginary parts. For this 
reason, it will be convenient for us to consider the following Hilbert operator: 



acting on functions in L P (1BL + ). Provided 1 < p < oo, the function Tiv + iv will then be in 
H p , with the presumption that v is an odd function. 

A natural type of problem in H 2 is that of approximating a given Hermitian function 
supported on two symmetric intervals. This means that one specifies a desired behavior 
in a certain frequency band and attempts to find, in some appropriate sense, an optimal 
approximation in H 2 . Without further restrictions, such a problem makes little sense, 
because it is easy to see that approximations can be made with arbitrary precision in 
L 2 norm. It may be reasonable to prescribe bounds for the norm of the approximating 
function; see for instance the work of M. G. Krem and P. Ya. Nudel'man [3] for 
interesting results along such lines. In the present note, we take a different route. We shall 
require the imaginary part of the function to be nonnegative for nonnegative arguments. 1 
We are interested in a specific problem of this general kind; it turns out to have an explicit 
and remarkably simple solution. 



We turn to the statement of the problem. For a finite interval I = [a, b] (0 < a < b) and 
every real number a we define the family of functions 



(Here and elsewhere we suppress the obvious "almost everywhere" provisions needed when 
considering pointwise restrictions.) We think of functions in K a (I), or more generally 
functions in L 2 (R + ), as the imaginary parts of Hermitian functions, and we view them 
therefore as odd functions on R. 

Our purpose is to give a parametrization of K a (I) and to solve the extremal problem 




2. Results 
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This reflects the passivity condition for our electromagnetic medium. 
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where xi denotes the characteristic function of /. We will show that the extremal problem 
has the following explicit solution: 

_ b 2 - a 2 
lab 

We note that the quantity on the right is invariant under dilations si = [sa,sb]. This is 
as it should be since v(t) is in K-\{I) if and only if v(t/s) is in K^i(sl) for s > 0. It will 
become clear that the extremal value A is not attained by any function in K_i(I). We will 
also see that the problem is insensitive to which LP norm we choose to minimize. 

Clearly, the corresponding extremal problem for K a (I) has solution |a|A when a < 0. 
However, if a > 0, the extremal problem is uninteresting and has solution 0. Thus the sign 
in the relation Tiv{t) = a for t 6 I matters in a decisive way. 2 

The following lemma is basic for our parametrization of K a (I). 

Lemma 1. A real-valued function in L 2 (R + \ I) is the restriction to R + \ / of at most one 
real-valued function v in L 2 (M. + ) such that T~Cv(t) is constant on I. 

Proof. We assume two real- valued functions V\ and v 2 in L 2 (R + ) coincide off / and are such 
that both U\(t) = Hvi(t) and u 2 (t) = 7iv 2 {t) are constant on I. If we set c = u 2 (t) — U\(t) 
for t in J, then the function h = [(ui — u 2 + i{v\ — v 2 ) + c] 2 will be real for real arguments. 
A change of variables argument shows that then h(i(l + z)/(l — z)) belongs to H l of the 
unit disk. But a function in H l can be real only if it is a constant. Clearly, h can be a 
constant only if u\ — u 2 + i{v\ — v 2 ) = 0. ■ 

We note that the assumption that v is in L? is essential for this lemma; the proof would 
break down if we assumed, say, that v belonged to some LP for p < 2. 
The following function will play an essential role in what follows: 



V 7 — i i 

a/ z 2 — b 2 y z 2 — a 2 

This function, which is taken to be positive for real arguments x > b, is analytic in the slit 
plane C \ ([—b, —a] U [a, b}). For real arguments a < \x\ < b we define a(x) by extending it 
continuously from the upper half-plane. Thus o~{x) takes values on the negative imaginary 
half-axis when x is in (a, b) and on the positive imaginary half-axis when — x is in (a, b), and 
otherwise it is real for real arguments. The key point, besides the symmetry a(— x) = cr(x), 
is that a provides a means for switching between real and imaginary when switching off 
and on /. 

The following is our main result. 



The case a < — 1 corresponds to the interesting physical phenomenon of negative refraction, which has 
received considerable attention in recent years. Artificial, negatively refracting materials, called metama- 
terials, have been realized in the microwave range 0, building on previous theoretical ideas As 
explained in [S] , the solution to our extremal problem provides a bound for the loss of negatively refracting 
materials when the real part of the refractive index is constant in a finite bandwidth. 
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Theorem. A nonnegative function v in L 2 (M. + ) is in K a (I) if and only if the following 
three conditions hold: 



(1) / v(t)\a(t)\dt < oo 

Jr+\i 



(2) - / tv(t)a(t)dt = a 



n Jr+\i 

(3) v(x) =H((1- X i)(tv) (x)/\a(x)\, x & I. 

Some remarks are in order before we give the proof of the theorem. 
The integrability condition is merely a slight growth condition at the endpoints of /; 
we may write it more succinctly as 



f a dt 
\ [v(a-t) +v(b + t)] —= < oo. 
Jo \t 



This condition ensures that the integral in (j2J and the Hilbert transform appearing in (jHJ) 
are both well-defined. 

At first sight, the theorem may not seem to give an explicit parametrization of K a (I). 
However, the Hilbert transform appearing in (jSJ) is given by 

If 2t 
H{{1- xi)av) (x) = - v(t)a(t)- -dt, 

7T J R +\j t — X 

and we observe that the integrand on the right is nonnegative whenever v(t) is nonnegative. 
Hence v(x) > for x off / implies v{x) > for x in /. This small miracle implies that 
K a (I) is parameterized by those nonnegative functions v in L 2 (R + \ J) for which (0) and 
(0) hold and such that 

\H((l -Xi)<™) {x)\ 2 \a(x)\~ 2 dx < oo. 



l 



By rephrasing this condition in more explicit terms (see Lemma 3 below), we arrive at the 
following corollary. 

Corollary. A nonnegative function v in L 2 (IR + \/) has an extension to a function in some 
class K a (I) if and only if the following condition holds: 

[u(a - t)u(a - r) + u(b + t)u(b + r)] ■ log ^ T ^ dtdr < oo. 



o J o 



tT 



The difference between (jj) and the condition above is the logarithmic factor, which 
means that the condition of the corollary is only a very slight strengthening of (JTJ). It is 
clear that for instance boundedness of v near the endpoints of I is more than enough. 
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We note that the integrand in (J2J) is negative to the left of I and positive to the right of 
/. This means that if a is negative, then 

2 r a 

\a\<- tv(t)\a(t)\dt, 



2 



7T _ _ 

with equality holding if v vanishes to the right of /. It follows that 

1 f a 2t \ct 
H{{l-xi)vv){x)>- / v{t)a{t)- 2 ^dt> 

7T / n X X X 



we may come as close as we wish to this lower bound by choosing any suitable v supported 
on a small set sufficiently close to 0. Hence our extremal problem (corresponding to 
a = — 1) has solution 



J(b 2 -x 2 )(x 2 -a 2 ) b 

A = max 



2 -a 2 



xei x 2 2ab ' 

as proclaimed above. We also observe that the same function 1/(x 2 |<t(x)|) would give the 
infimum for the L p norm over I for any other value of p > 0. 
If, on the other hand, a is positive, we have instead 

2 f°° 

a < - / tv(t)\a(t)\dt, 
n Jb 

with equality holding if v vanishes to the left of J. In this case, arguing in the same fashion 
as above, we find that we can get Tt ((1 — Xi) av ) ( x ) as small as we please by letting v be 
supported on a set sufficiently far to the right of /. 

3. Proofs 

We now turn to the proof of the theorem and its corollary. We will rely on Lemma 1 
and two additional lemmas. 

Lemma 2. For every t in (0, a) U (b, oo) the function 

AM = ;Aj (i - Tl) - 2ta{t) 

t 2 — x 2 \ a(x) J 
is in H p for p > 1/2, and the following estimates hold 

WMU < — =Jq=== t<2b, 

\ft{x)\ 2 dt< ( ' 2 



2b 



X\ + l' 



where the constants C% and C2 only depend on a and b. 

Proof. It is immediate that f t belongs to H p for p > 1/2 because the isolated singularities 
±t are removable and ft{z) = 0(z~ 2 ) when z — > 00. The norm estimates follow from 
elementary calculations. ■ 



6 KRISTIAN SEIP AND JOHANNES SKAAR 

Lemma 3. A nonnegative function v in L 2 (W + \ I) satisfies 

\H({1 - Xi)o-v)(x)\ 2 \o-(x)\~ 2 dx < oo 
if and only if the following condition holds: 

[v(a - t)u{a - r) + u(b + t)v(b + r)] I lo s( t + r ) I dtdr < ^ 



o Jo 



tT 



Proof. The necessary and sufficient condition for square-integrability at the left end-point 
of / is that 

via — t) \ 2 

-dt x dx < oo. 



o jo 



'o Vi(t + x) 

By Fubini's theorem, we may interchange the order of integration so that this condition 
becomes 

' . . , |log(t + r)| , , 
v[a — t)v[a — t) = dtdr < oo. 

\JtT 

Combining this with the corresponding condition at the right end-point of /, we arrive at 
the condition of the lemma. ■ 

The theorem is now proved in the following way. We assume first that we are given a 
nonnegative function v in L 2 (R + ) satisfying (JIJ, (J2J), and (jBJ)- We claim that the function 

fi(x) = - f v(t)f t (x)dt 

n J (0,a)U(b,2b) 

is in H . Indeed, by Lemma 2 and Fubini's theorem, 

«(*) 



\\h\V<^f 

n Ja 



, ^=dt, 

(0,o)U(6,26) y/\(a — t)(b — t)\ 

and the integral on the right is bounded thanks to On the other hand, 

f 2 ( x ) = I f v(t)f t {x)dt 



n J2b 

is in H p for p > 2, because by the Cauchy-Schwarz inequality and Lemma 2 we have 

( - / 2 " ' 1 



\f 2 (x)\ p dx < -^-llull? / -r,dx. 

The imaginary part of f\ + fa is supported by / and equals — v there, in view of (jHJ). Since 
v is assumed to be in L 2 , it follows that fi + f 2 is in fact in H 2 . 

We set / = 7i((l — Xi) v ) + *(1 ~ Xi) v -> which is a function in H 2 . We observe that the 
imaginary part of / — (/i + f 2 ) equals v and that its real part equals a on J, when taking 
into account (J2J). So we have proved that the given v is indeed in K a (I). 

We now prove the necessity of the three conditions of the theorem. So assume we 
are given a nonnegative function v in L 2 (M + ) belonging to some class K a {I). Setting 
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I E — [a + e, b — e], we see that v also belongs to K a (I £ ) whenever < e < {b — o)/2. Bnt 
since the Hilbert transform of v is constant near the endpoints of I £ , it follows that 



/ 



v(t)\a £ (t)\dt < oo, 

•/R+\I 6 

where now 

1 



y/t*-(b-e) 2 y/t 2 -(a + e) 2 ' 
We claim that this means that 

(4) v(x)=H((l- X i a )<rev)(x)/\a e (x)\, 

provided x is in (a + e, b — e). Indeed, by Lemma 1, it is enough to verify that 

(1 - Xle(x))v(x) + XleWH ((1 - Xle)°sV) ( x) / \ G £ (x) \ 

is in K a (I e ) for some a. Since, in view of Lemma 3, the function on the right-hand side of 
(j3J) is square-integrable on I £ , the claim follows by repeating the argument in the first part 
of the proof. 

We may view TC ((1 — Xie) a e v ) ( x ) as the L 1 norm of the function 

1 2t 

h x>£ (t) = - ((1 - X i e (t))a e (t)v(t)) 

Then (jlj) says that ||/i x ,e||i —* v(x)\a{x)\ when e —>■ 0. Since we also have that 

1 2t 
h x ,M ->-((!- Xi(t))a(t)v(t)) 



71 t 2 — X 2 

for every t, we obtain 

v(x)=H((l-xi)™)(x)/\a(x)\ 
for every x in (a, 6). By a similar argument, we find that 

a = lim- / tv(t)a E (t)dt = - I tv(t)a(t)dt. 

The necessity of (jlj has already been observed; without it we would reach the contradictory 
conclusion that v(x) = oo for almost every x G (a, 6). 

We finally note that the corollary is an immediate consequence of the theorem and 
lemmas 1 and 3. 
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